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ABSTRACT 

Many high-state non-magnetic cataclysmic variables (CVs) exhibit blue-shifted ab¬ 
sorption or P-Cygni profiles associated with ultraviolet (UV) resonance lines. These 
features imply the existence of powerful accretion disk winds in CVs. Here, we use our 
Monte Carlo ionization and radiative transfer code to investigate whether disk wind 
models that produce realistic UV line profiles are also likely to generate observation- 
ally significant recombination line and continuum emission in the optical waveband. 
We also test whether outflows may be responsible for the single-peaked emission line 
prohles often seen in high-state CVs and for the weakness of the Balmer absorption 
edge (relative to simple models of optically thick accretion disks). We find that a stan¬ 
dard disk wind model that is successful in reproducing the UV spectra of CVs also 
leaves a noticeable imprint on the optical spectrum, particularly for systems viewed at 
high inclination. The strongest optical wind-formed recombination lines are Ha and 
He II A4686. We demonstrate that a higher-density outflow model produces all the 
expected H and He lines and produces a recombination continuum that can hll in the 
Balmer jump at high inclinations. This model displays reasonable verisimilitude with 
the optical spectrum of RW Trianguli. No single-peaked emission is seen, although we 
observe a narrowing of the double-peaked emission lines from the base of the wind. 
Finally, we show that even denser models can produce a single-peaked Ha line. On the 
basis of our results, we suggest that winds can modify, and perhaps even dominate, 
the line and continuum emission from CVs. 

Key words: (stars:) novae, cataclysmic variables - accretion, accretion discs - stars: 
winds, outflows - line: profiles - radiative transfer - methods: numerical 


1 INTRODUCTION 

Cataclysmic variables (CVs) are systems in which a white 
dwarf accretes matter from a donor star via Roche-lobe 
overflow. In non-magnetic systems this accretion is medi¬ 
ated by a Keplerian disk around the white dwarf (WD). 
Nova-like variables (NLs) are a subclass of CVs in which 
the disk is always in a relatively high-accretion-rate state 
{M 10“® Mq yr“^). This makes NLs an excellent labo¬ 
ratory for studying the properties of steady-state accretion 
disks. 

It has been known for a long time that winds ema¬ 
nating from the accretion disk are important in shaping 
the u l traviolet (UV) spectra o f high-state CVs (iHeap et al.l 
Il978l : iGreenstein fc oSll982l l. The most spectacular evi¬ 
dence for such outflows are the P-Cygni-like profiles seen 
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in UV resonance lines such as C iv A1550 (see e.g. Cor¬ 
dova & Mason 1982). Considerable effort has been spent 
over the years on understanding and modelling these UV 
features (e.g. Drew & Verbunt 1985; Mauche & Raymond 
1987; Drew 1987; Shlosman & Vitello 1993; [hereafter SV93]; 
Knigge, Woods & Drew 1995; Knigge & Drew 1997; Knigge 
et al. 1997; Long & Knigge 2002 [hereafter LK02[, Noebauer 
et al. 2010; Puebla et al. 2011). The basic picture emerging 
from these efforts is of a slowly accelerating, moderately 
collimated bipolar outflow that carries away ~ 1% — 10% of 
the accreting material. State-of-the-art simulations of line 
formation in this type of disk wind can produce UV line 
profiles that are remarkably similar to observations. 

Much less is known about the effect of these outflows 
on the optical spectra of high-state CVs. These spectra are 
typically characterized by H and He emission lines super¬ 
posed on a blue continuum. In ma ny cases, and particu- 
larly in the SW Sex subclass of NLs jHonevcutt et aHll986l : 
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Figure 1. Optical spectra of three nova-like variables: RW Sex (top; Beuermann et al. 1992), IX Vel (top middle; A. F. Pala & B. 
T. Gaensicke, private communication) and RW Tri in and out of eclipse (bottom two panels; Groot et al. 2004). The data for RW Sex 
and RW Tri were digitized from the respective publications, and the IX Vel spectrum was obtained using the XSHOOTER spectrograph 
on the Very Large Telescope on 2014 October 10. These systems have approximate inclinations of 30°, 60° and 75° (see section 5.4) 
respectively. The trend of increasing Balmer line emission with inclination can be seen. In RW Tri strong single-peaked emission in the 
Balmer lines is seen even in eclipse, indicating that the lines may be formed in a spatially extensive disk wind, and there is even a 
suggestion of a (potentially wind-formed) recombination continuum in the eclipsed spectrum. We have attempted to show each spectrum 
over a similar dynamic range. 


IPhillon fc R.utt^ll99!Tl ). these lines are single-peaked. This 
is contrary to theoretical expectations for lines formed in 
accretion di s ks, which are predict ed to be double-peaked 
llSmakifloSlI : iHorne &: Marshlll986l) . Low-state CVs (dwarf 
novae in quiescence) do, in fact, exhibit such double-peaked 
lines (iMarsh &: Hornelll990l ). 

Murray & Chiang (1996, 1997; hereafter referred to col¬ 
lectively as MC96) have shown that the presence of disk 
winds may offer a natural explanation for the single-peaked 
optical emission lines in high-state CVs, since they can 
strongly affect the radiative transfer of line photons. Strong 
support for a significant wind contribution to the optical 
emission lines comes from observations of eclipsing systems. 
There, the single-peaked lines are often only weakly eclipsed, 
and a significant fraction of the line flux remains visible 
even near mid-eclipse (e.g. lBaotista et aI.ll200Cl| : ICroot et al.l 
I 2 OO 4 I ). This points to line formation in a spatially extended 
region, such as a disk wind (see Fig. [T|). Further evidence 
for a wind contribution to the optical lines comes from iso¬ 
lated observations of P-Cygni-like lin e profiles even in opti¬ 
cal lines, such as Hq and He l A58 76 (IPattersqn et al.| ll99B : 
iRinewald fc Navlorj|l998l : iKafka fc Honevcutd 20041 ). 

Could disk winds also have an impact on the 
UV/optical continuum of high-state CVst This contin¬ 
uum is usually thought to be dominated by the accre¬ 
tion disk and modelled by splitting the disk into a set 
of concentric, optically thick, non-interacting annuli fol¬ 
lowing the stand ard Teff(R) oc r adial tempera¬ 

ture distribution llShakura fc Sunvaevlll973 ). In such mod¬ 
els, each annulus is taken to emit either as a black- 
body or, perhaps more realistically, as a stellar/disk at¬ 


mosphere model llSchwarzenberg-Czernv fc R6zvczkalll977l : 
I Wadelfr984l . 1 198§ ) . In the latter case, the local surface grav¬ 
ity, log g{R), is assumed to be set solely by the accreting 
WD, since self-gravity is negligible in CV disks. 

Attempts to fit the observed spectral energy distribu¬ 
tions (SEDs) of high-state CVs with such models have met 
with mixed success. In particular, the SEDs predicted by 
most stellar/disk atmosphere models a re too blue in the UV 
dWadd 1198^ iLong et aH Il99ll . 1 1994 iKnigge et al.l Il998ah 


and exhibit strong 

er-than-observed Balmer lumps in absorp- 

tion ( 

Waddll984 

Hauglll987l: iLa Doudll9894 iKnigge et al.l 

Il998a 

). One possible explanation for these problems is that 


these models fail to capture all of the relevant physics. In¬ 
deed, it has been argued that a self-consistent treatment can 
produce better agreement with observational data (e.g. Sha- 
viv et al. 1991; but see also Idan et al. 2010). However, an 
alternative explanation, suggested by Knigge et al. (1998b; 
see also Hassall et al. 1985), is that recombination contin¬ 
uum emission from the base of the disk wind might fill in the 
disk’s Balmer absorption edge and flatten the UV spectrum. 

Here, we carry out Monte Carlo radiative transfer sim¬ 
ulations in order to assess the likely impact of accretion disk 
winds on the optical spectra of high-state CVs. More specif¬ 
ically, our goal is to test whether disk winds of the type 
developed to account for the UV resonance lines would also 
naturally produce significant amounts of optical line and/or 
continuum emission. In order to achieve this, we have im¬ 
plemented the ‘macro-atom’ approach developed by Lucy 
(2002, 2003) into the Monte Carlo ionization and radiative 
transfer code described by LK02 (a process initiated by Sim 
et al. 2005; hereafter SDL05). With this upgrade, the code is 
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able to deal correctly with processes involving excited levels, 
such as the recombination emission produced by CV winds. 

The remainder of this paper is organized as follows. In 
Section 2, we briefly describe the code and and the newly im¬ 
plemented macro-atom approach. In Section 3, we describe 
the kinematics and geometry of our disk wind model. In Sec¬ 
tion 4, we present spectra simulated from the benchmark 
model employed by LK02, and, in Section 5, we present a 
revised model optimized for the optical waveband. In Sec¬ 
tion 6, we summarize our findings. 


2 PYTHON: A MONTE CARLO IONIZATION 
AND RADIATIVE TRANSFER CODE 

Python is a Monte Carlo ionization and radiative trans¬ 
fer code whi ch uses the Sobolev approximation to treat lin e 
transfer fe.g. lSobolevlll957l . ll960l : [Rvbicki &: Hummeilll978li . 
The code has already been described extensively by LK02, 
SDL05 and Higginbottom et al. (2013; hereafter H13), so 
here we provide only a brief summary of its operation, fo¬ 
cusing particularly on new aspects of our implementation of 
macro-atoms into the code. 


2.1 Basics 

Python operates in two distinct stages. First, the ionization 
state, level populations and temperature structure are calcu¬ 
lated. This is done iteratively, by propagating several popu¬ 
lations of Monte Carlo energy quanta (‘photons’) through a 
model wind. The geometric and kinemetic properties of the 
outflow are specified on a pre-defined spatial grid. In each of 
these iterations (‘ionization cycles’), the code records esti¬ 
mators that characterize the radiation field in each grid cell. 
At the end of each ionization cycle, a new electron temper¬ 
ature is calculated that more closely balances heating and 
cooling in the plasma. The radiative estimators and updated 
electron temperature are then used to revise the ionization 
state of the wind, and a new ionization cycle is started. The 
process is repeated until heating and cooling are balanced 
throughout the wind. 

This converged model is then used as the basis for the 
second set of iterations (‘spectral cycles’). In these, the emer¬ 
gent spectrum over the desired spectral range is synthesized 
by tracking populations of energy packets through the wind 
and computing the emergent spectra at a number of user- 
specified viewing angles. 

Python is designed to operate in a number of differ¬ 
ent regimes, both in terms of the scale of the system and in 
terms of the characteristics of the underlying radiation field. 
It was originally developed by LK02 in order to model the 
UV spectra of CVs with a simple biconical disk wind model. 
SDL05 used the code to model Brackett and Pfund line pro¬ 
files of H in young-stellar objects (YSOs). As part of this ef¬ 
fort, they implemented a ‘macro-atom’ mode (see below) in 
order to correctly treat H recombination lines with Python. 
Finally, H13 used Python to model broad absorption line 
(BAL) QSOs. For this application, an improved treatment 
of ionization was implemented, so that the code is now capa¬ 
ble of dealing with arbitrary photo-ionizing SEDs, including 
non-thermal and multi-component ones. 


2.2 Ionization and Excitation: ‘Simple Atoms’ 

Prior to SDL05, the relative ionization fractions for all 
atomic species were estimated via the modified Saha equa¬ 
tion (Mazzali & Lucy 1993) 

= + ( 1 ) 

Here, the ‘starred’ term on the right represents abundances 
computed with the Saha equation at temperature Tr, but 
using partition functions from the dilute blackbody approx¬ 
imation. W is an effective dilution factor, ^ is the fraction of 
recombinations going directly to the ground state, and Tr 
and Te are the radiation and electron temperatures, respec¬ 
tively. This simple ionization scheme produces reasonable 
results when the photoionizing SED can be approximated 
by a dilute blackbody. This is the case for high-state CVs. 
(As noted above, an improved, but more complex treatment 
of ionization that is appropriate for more complex SEDs is 
described in H13.) 

Similarly, the relative excitation fractions within each 
ionization stage of a given species were estimated via a mod¬ 
ified (dilute) Boltzmann equation, 

where rijk is the population of level k in ionic stage j, Ek is 
the energy difference between level k and the ground state, 
Qk is the statistical weight of level k and zATr) is the par¬ 
tition function of ionic stage j. 

Finally, Python originally modelled all bound-bound 
processes as t ransitions within a simple two-level atom (e.g. 
lMihalaslll982h . This framework was used for the treatment 
of line transfer and also for the line heating and cooling cal¬ 
culations (see LK02). The approximation works reasonably 
well for resonance lines, such as C iv A1550, in which the 
lower level is the ground state. However, it is a poor ap¬ 
proximation for many other transitions, particularly those 
where the upper level is primarily populated from above. 
Thus an improved method for estimating excited level pop¬ 
ulations and simulating line transfer is needed in order to 
model recombination lines and continua. 


2.3 Ionization and Excitation: Macro-Atoms 


Lucy (2002, 2003; hereafter L02, L03) has shown that it is 
possible to calculate the emissivity of a gas in statistical 
equilibrium accurately by quantising matter into ‘macro¬ 
atoms’, and radiant and kinetic energy into indivisible en¬ 
ergy packets (r- and k- packets, respectively). His macro¬ 
atom scheme allows for all possible transition paths from 
a given level and provides a full non-local thermodynamic 
equilibrium (NLTE) solution for the level populations based 
on Monte Carlo estimators. The macro-atom techniqu e has 
alrea dy been used to model Wolf-Rayet star winds (jSiml 
l2004h , AGN disk winds (lSjm et_^J|2008^;|Thtmnl^ 201^, 
supernovae dKromer fc Siml l2009l : iKerzendorf fc SimI 1201411 
and YSOs (SDL05). A full description of the approach can 
be found in L02 and L03. 

Briefly, macro-atom NLTE level populations and ioniza¬ 
tion fractions are calculated by solving the statistical equi- 
Ibrium equations between each pair of levels. In the frame¬ 
work of the Sobolev escape probability formalism (Rybicki & 
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Hummer 1978; L02; Sim 2004), the bound-bound excitation 
rate, TZiu, in an ion is given by 

(3) 

where u and I denote the upper and lower levels, C repre¬ 
sents the collisional rate coefficients, and B is the usual Ein¬ 
stein coefficient. J^st is the Monte Carlo estimator for the 
mean intensity impinging on the Sobolev region, weighted by 
an angle-dependent escape probability, given by (ISimll2004l l 

Jest — ^ / 'Wi / J \ ’ (‘^) 

Atvi/qV ^ Ts^i {dv/ds)i 

Here w is the photon weight (in luminosity units), is 
the line frequency, dv/ds is the velocity gradient and Ts is 
the Sobolev optical depth. The sum is over all photons that 
come into resonance with the line, and thus represents an 
integral over solid angle. The corresponding de-excitation 
rate is then 


t^ul — Bul^uJest 4 - CulT^u^e^ ( 5 ) 


where A is the usual Einstein coefficient. The quantity Piu is 
the angle-averaged probability that a given line photon will 
escape the Sobolev region. 

In our implementation of the macro-atom approach, we 
also explicitly take into account the photoionization and col¬ 
lisional ionization rates between a lower level, I, and the con¬ 
tinuum (or, in the case of ions with more than one bound 
electron, the ground state of the upper ion), k. 


B-Ik, — Til 



hu 


dv -\- Ci^nin^. 


( 6 ) 


Here, is the photoionization cross section, and Ju is the 
mean intensity. The corresponding recombination rate is 
given by 


'B-tzl — TTlfzlTli^Tle 4 - CtzlTlt^Tle^ 


(7) 


where UkI is the radiative recombination coefficient to level 
1. This treatment means that radiative and collisional rates 
to and from all levels are considered when calculating both 
the ionization state and the level populations, although we 
negl ect ionization dir e ctly t o excited levels of the upper ion. 
The Ivan Regemorteii (Il962l ~) approximation is used for col¬ 
lisional transitions. This means that collisions between ra- 
diatively forbidden transitions are not taken into account 
when one splits levels ito I- and s-subshells, as well as prin¬ 
cipal quantum number, n (as we have done with He l; see 
section 12.511 . Although this approximation is, in general, a 
poor one, the effect is second order in the physical regime 
where recombination lines are formed in our models. This is 
because bound-free processes are dominant in determining 
level populations and emissivities. We have verified that this 
is indeed the case in the He l emission regions in our models. 


2.4 Ionization and Excitation: A Hybrid 
Approach 

SDL05 implemented a macro-atom treatment of H in 
Python and used this to predict the observable proper¬ 
ties of a pure H wind model for YSOs. Our goal here is to 
simultaneously model the optical and ultraviolet spectra of 
high-state CVs. Since the optical spectra are dominated by 


H and He recombination lines, both of these species need to 
be treated as macro-atoms. The UV spectra, on the other 
hand, are dominated by resonance lines associated with met¬ 
als. This means we need to include these species in our mod¬ 
els, but they can be treated with our (much faster) simple- 
atom approach. We have therefore implemented a hybrid 
ionization and excitation scheme into Python. Any atomic 
species can now be treated either in our simple-atom ap¬ 
proximation or with the full macro-atom machinery. In our 
CV models, we treat H and He as macro-atoms and all met¬ 
als as simple-atoms. Species treated with either method are 
full taken into account as sources of both bound-free opacity 
and line opacity, and contribute to the heating and cooling 
balance of the plasma. 


2.5 Atomic Data 

We generally use the same atomic data as H13, which is 
an updated version of that described by LK02. In addition, 
we follow SDL05 in treating H as a 20-level atom, where 
each level is defined by the principal quantum number, n. 
For the macro-atom treatment of He, we have added the ad- 
ditional level and li ne information required from Topbase 
llBadnell et akll^OOSh . He II is treated in much the same way 
as H, but with 10 levels. He l has larger energy differentials 
between different 1-subshells and triplet and singlet states. 
Thus, we still include levels up to n = 10, but explicitly 
treat the I and s sub-orbitals as distinct levels instead of as¬ 
suming they are perfectly ‘1-mixed’. This allows us to model 
the singlet and triplet He l lines that are ubiquitous in the 
optical spectra of CVs (e.g. iDhillonlflo^ . 


2.6 Code Validation and Testing 

Python has been tested against a number of radiative 
transfer and photoionization codes. LK02 and H13 con- 
ducted comparisons of ionization balance with Cloudy 
dFerland et al.l l2013ll . demonstrating excellent agreement. 
We have also carried out comparisons of ionization and 
spectral synthe sis with the supernova cod e Tardis. Tardis 
is described by iKerzendorf fc SimI (l2014h . and the spectral 
comparisons can be found therein. For the effort reported 
here, we have additionally carried out tests of the macro¬ 
atom scheme in Python. Fig.[2]shows two of these tests. In 
the top panel, we compare the Balmer series emissivities as 
predicted by Python in the 1-mixed Case B limit against 
the analytical calculations bv ISeatonI (Il959l l. In the bottom 
panel, we compare Python and Tardis predictions of He 
I level populations for a particular test case. Agreement is 
excellent for both H and He. 


3 DESCRIBING THE SYSTEM AND ITS 
OUTFLOW 

Python includes several different kinematic models of ac¬ 
cretion disk winds, as well as different options for describ¬ 
ing the physical and radiative properties of the wind-driving 
system under consideration. Most of these features have al¬ 
ready been discussed by LK02 and H13, so below we only 
briefly recount the key aspects of the particular system and 
wind model used in the present study. 
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Figure 2. Top Panel: ‘Case B’ Balmer decrements computed 
with Python compared to analytic calculations by Seaton (1959). 
Both calculations are calculated at Tg = 10,000K. (see Oster- 
brock 1989 for a discussion of this commonly used approxima¬ 
tion). Bottom Panel: a comparison of He I level populations (the 
most complex ion we currently treat as a macro-atom) between 
Python and Tardis models. The calculation is conducted with 
physical parameters of rig = 5.96 x lO'^ cm~^, Tg = 30, 600K, 
Tji = 43, 482K and W = 9.65 x 10“^. Considering the two codes 
use different atomic data and Tardis, unlike Python, currently 
has a complete treatment of collisions between radiatively forbid¬ 
den transitions, the factor of < 2 agreement is encouraging. 


3.1 Wind Geometry and Kinematics 

We adopt the kinematic disk wind model developed by 
SV93. A schematic of this model is shown in Fig. [3] In this 
parametrization, a smooth, biconical disk wind emanates 
from the accretion disk between radii rmin and Vmax^ The 
covering fraction of the outflow is also controlled by the in¬ 
ner and outer opening angles of the wind, Omin and Omax, 
and the launch angle of the other streamlines is given by 

e{ro) = + ( 0 ™. - , ( 8 ) 

\'Omax fmin / 

where ro is the launch radius of the streamline. 

The poloidal (non-rotational) velocity field of the wind, 
vi, is given by 

, r . ^ , d/R^r 

VI = V0 + boo(ro) - uo] (9) 

where I is the poloidal distance along a particular wind 
streamline. The terminal velocity along a streamline, Uoo, 
is set to a fixed multiple of Uesc, the escape velocity at the 
launch point. The launch velocity from the disk surface, vo, 
is assumed to be constant (set to 6 km s“^). Once the wind 
is launched, it accelerates, reaching half of its terminal ve¬ 
locity at I — Rv The velocity law exponent a controls how 
quickly the wind accelerates. Larger values of a cause the 
main region of acceleration to occur close to Rv, whereas 
smaller values correspond to fast acceleration close to the 
disk (see Fig. 0. The rotational velocity is Keplerian at 
the base of the streamline and we assume conservation of 
specific angular momentum, such that 

( 10 ) 



Figure 3. Cartoon illustrating the geometry and kinematics of 
the benchmark CV wind model. 



Figure 4. The adopted poloidal velocity law for various values 
of the acceleration exponent, a. 


where Vk = (GMwn/ro)^^^- 

The density at position (r, z) in the wind, p(r,z), is 
calculated from the mass continuity equation, yielding 


r dr Vz{r,2) 


( 11 ) 


Here, Vz is the vertical velocity component and, following 
SV93, <(>(ro) is the local mass-loss rate per unit area at ro, 
defined as 


(f)(ro) oc MniindTo cos[0(ro)]. (12) 

We adopt A = 0 and normalize 4>{ro) by matching its integral 
over both sides of the disk to the user-specified total mass- 
loss rate, M^ind- 


3.2 Sources and Sinks of Radiation 

The net photon sources in our CV model are the accretion 
disk, the WD and, in principle, a boundary layer with user- 
defined temperature and luminosity. All of these radiating 
bodies are taken to be optically thick, and photons strik¬ 
ing them are assumed to be destroyed instantaneously. The 


v^r = VkVo, 
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secondary star is not included as a radiation source, but 
is included as an occulting body. This allows us to model 
eclipses. Finally, emission from the wind itself is also ac¬ 
counted for, but note that we assume the outflow to be in 
radiative equilibrium. Thus all of the heating of the wind, 
as well as its emission, is ultimately powered by the radia¬ 
tion field of the net photon sources in the simulation. In the 
following sections, we will describe our treatment of these 
system components in slightly more detail. 


3.2.1 Accretion Disk 


Python has some flexibility when treating the accretion 
disk as a source of photons. The disk is broken down into 
annuli such that each annulus contributes an equal amount 
to the bolometric luminosity. We take the disk to be geomet¬ 
rically thin, but optically thi ck, and thus adopt the tem per- 
ature profile of a standard IShakura fc SunvaevI (Il973f) a- 
disk. An annulus can then be treated either as a blackbody 
with the corresponding effective temperature or as a stel¬ 
lar atmosphere model with the appropriate surface gravity 
and effective temperature. Here, we use blackbodies dur¬ 
ing the ionization cycles and to compute our Monte Carlo 
estimators. However, during the spectral synthesis stage of 
the simulation we use stellar atmosphere models. This pro¬ 
duces more realistic model spectra and allows us to test if 
recombination emission from the wind base can fill in the 
Balmer jump, which is always in absorption in these models. 
Our synthetic stellar atmospher e spectra are c alculated with 
SynspecQ from either Kurucz llKurucjll991) atmospheres 
(for T eff ^ 50,000 K) or from TLUSTY ( Hubenv fc Lanj 
Il995^ models (for T^ff > 50,000 K). 


3.2.2 White Dwarf 

The WD at the center of the disk is always present as a 
spherical occulting body with radius Rwd in Python CV 
models, but it can also be included as a source of radia¬ 
tion. In the models presented here, we treat the WD as a 
blackbody radiator with temperature Two and luminosity 
Lwd = 47ri?^y £) £). 


3.2.3 Boundary Layer 

It is possible to include radiation from a boundary layer 
(BL) between the disk and the WD. In Python, the BL 
is described as a blackbody with a user-specified effective 
temperature and luminosity. In the models presented here, 
we have followed LK02 in setting the BL luminosity to zero. 
However, we have confirmed that the addition of an isotropic 
BL with Lbl = 0.5Lacc and temperatures in the range 
80 kK ^ Tbl ^ 200 kK would not change any of our main 
conclusions. 


3 . 2.4 Secondary Star 

The donor star is included in the system as a pure radiation 
sink, i.e. it does not emit photons, but absorbs any pho¬ 
tons that strike its surface. The secondary is assumed to be 

^ http://nova.astro.umd.edu/Synspec43/synspec.html 


Model Parameters 


Parameter 

Model A 

Model B 

Mwd 

0.8 Mq 


Rwd 

7 X 10® cm 


Twd 

40, 000 K 


M2 

- 

0.6 Mq 

9 

- 

0.75 

^orb 

- 

5.57 hr 

a 

- 

194.4 Rwd 

R2 

- 

69.0 Rwd 

^acc 

10~® MQyr~^ 


^wind 

10“® MQyr~^ 


^ min 

4 Rwd 


^ max 

12 Rwd 


Tdisfc (max) 

34.3 Rwd 


^min 

20.0° 


^max 

65.0° 


7 

1 



3 Ve.sc 


Rv 

100 Rwd 

142.9 Rwd 

a 

1.5 

4 


Table 1. Parameters used for the geometry and kinematics of the 
benchmark CV model (model A), which is optimized for the UV 
band, and a model which is optimized for the optical band and 
described in section 5 (model B). For model B, only parameters 
which are altered are given - otherwise the model A parameter is 
used. Parb is the orbital period (the value for RW Tri from Walker 
1963 is adopted, see section 5.4) and R 2 is the radius of a sphere 
with the volume of the secondary’s Roche lobe. Other quantities 
are defined in the text or Fig. [3] Secondary star parameters are 
only quoted for model B as we do not show eclipses with the 
benchmark model (see section 5.4). 


Roche-lobe filling, so its shape and relative size are defined 
by setting the mass ratio of the system, q = M 2 /Mwd. The 
inclusion of the donor star as an occulting body allows us to 
model eclipses of the disk and the wind. For this purpose, we 
assume a circular orbit with a semi-major axis a and specify 
orbital phase such that ‘Lorb = 0 is the inferior conjunction 
of the secondary (i.e. mid-eclipse for i ~ 90°). 


4 A BENCHMARK DISK WIND MODEL 

Our main goal is to test whether the type of disk wind model 
that has been successful in explaining the UV spectra of CVs 
could also have a significant impact on the optical continuum 
and emission line spectra of these systems. In order to set 
a benchmark, we therefore begin by investigating one of the 
fiducial CV wind models that was used by SV93 and LK02 
to simulate the UV spectrum of a typical high-state system. 
The specific parameters for this model (model A) are listed 
in Table 1. A key point is that the wind mass-loss rate in 
this model is set to 10% of the accretion rate through the 
disk. The inner edge of the wind (r^in) is set to 4 Rwd 
following SV93. The sensitivity to some of these parameters 
is briefly discussed in section 5. 
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Figure 5. The physical properties of the wind - note the logarithmic scale. Near the disk plane the wind is dense, with low poloidal 
velocities. As the wind accelerates it becomes less dense and more highly ionized. The dominant He ion is almost always He HI, apart 
from in a small portion of the wind at the base, which is partially shielded from the inner disk. 


4.1 Physical Structure and Ionization State 

Fig. [5] shows the physical and ionization structure of the 
benchmark disk wind model. The ionization parameter 
shown in the bottom right panel is given by 


Itt /■“ J^du 

nnC 713,6eV 


(13) 


where uh is the local number density of H, and v denotes 
photon frequency. The ionization parameter is a useful mea¬ 
sure of the ionization state of a plasma, as it evaluates the 
ratio of the number density of ionizing photons to the local 
H density. 

There is an obvious drop-off in density and temperature 
with distance away from the disk, so any line formation pro¬ 
cess that scales as - i.e. recombination and collisionally 
excited emission - should be expected to operate primar¬ 
ily in the dense base of the outflow. Moreover, a compari¬ 
son of the rotational and poloidal velocity fields shows that 


rotation dominates in the near-disk regime, while outflow 
dominates further out in the wind. 

The ionization equation used in the ‘simple atom’ ap¬ 
proach used by LK02 (see section should be a reason¬ 
able approximation to the photoionization equilibrium in 
the benchmark wind model. Even though the macro-atom 
treatment of H and He does affect the computation of the 
overall ionization equilibrium, we would expect the resulting 
ionization state of the wind to be quite similar to that found 
by LK02. The bottom panels in Fig. [5] confirm that this is 
the case. In particular. He is fully ionized throughout most 
of the outflow, except for a small region near the base of the 
wind, which is shielded from the photons produced by the 
hot inner disk. In line with the results of LK02, we also find 
that Civ is the dominant C ion throughout the wind, result¬ 
ing in a substantial absorbing column across a large range 
of velocities. As we shall see, this produces the broad, deep 
and blue-shifted Civ A1550 absorption line that is usually 
the most prominent wind-formed feature in the UV spectra 
of low-inclination nova-like CVs. 
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Figure 6 . UV (left) and optical (right) synthetic spectra for model A, our benchmark model, computed at sightlines of 10, 27.5, 45, 
62.5 and 80 degrees. The inset plots show zoomed-in line profiles for He ii A1640 and Ho. Double-peaked line emission can be seen in 
He II A1640, He ii A4686, Ha and some He I lines, but the line emission is not always sufficient to overcome the absorption cores from 
the stellar atmosphere models. The model also produces a prominent He ii A3202 line at high inclinations. 
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4.2 Synthetic Spectra 

We begin by verifying that the benchmark model still pro¬ 
duces UV spectra that resemble those observed in CVs. We 
do expect this to be the case, since the ionization state 
of the wind has not changed significantly from that com¬ 
puted by LK02 (see section ICT . The left column of panels 
in Fig.[6]shows that this expectation is met: all of the strong 
metal resonance lines - notably N v A1240, Si iv A1400 and 
C IV A1550 ” are present and exhibit clear P-Cygni pro¬ 
files at intermediate inclinations. In addition, however, we 
now also find that the wind produces significant Lya and 
He II A1640 emission lines. 

Fig. [6] (right-hand panel) and Fig. [7] show the corre¬ 
sponding optical spectra produced for the benchmark model, 
and these do exhibit some emission lines associated with H 
and He. We see a general trend from absorption lines to 
emission lines with increasing inclination, as one might ex¬ 
pect from our wind geometry. This trend is consistent with 
observations, as can be seen in Fig. 1. However, it is clear 
that this particular model does not produce all of the lines 
seen in observations of high-state CVs. The higher-order 
Balmer series lines are too weak to overcome the intrinsic 
absorption from the disk atmosphere, and the wind fails 
to produce any observable emission at low and intermedi¬ 
ate inclinations. This contrasts with the fact that emission 
lines are seen in the opt ical spectra of (for example) V3885 
Sgr (Hartley et al.l[20n3 ~) and IX Vel JBeuermann fc ThomasI 
ll99C . see also Fig. 1). 

The emissivity of these recombination features scales as 
p , meaning that they form almost entirely in the dense base 
of the wind, just above the accretion disk. Here, the velocity 
field of the wind is still dominated by rotation, rather than 
outflow, which accounts for the double-peaked shape of the 
lines. In principle, lines formed in this region can still be sin¬ 
gle peaked, since the existence of a poloidal velocity gradient 
changes the local escape probabilities (MC96). However, as 
discussed further in section 5.3, the radial velocity shear in 
our models is not high enough for this radiative transfer 
effect to dominate the line shapes. 

The Balmer jump is in absorption at all inclinations 
for our benchmark model. This is due to the stellar atmo¬ 
spheres we have used to model the disk spectrum; it is not 
a result of photoabsorption in the wind. In fact, the wind 
spectrum exhibits the Balmer jump in emission, but this is 
not strong enough to overcome the intrinsic absorption edge 
in the disk spectrum. This is illustrated in Fig. [8l which 
shows the angle-integrated spectrum of the system, i.e. the 
spectrum formed by all escaping photons, separated into the 
disk and wind contributions. Even though the wind-formed 
Balmer recombination continuum does not completely fill in 
the Balmer absorption edge in this model, it does already 
contribute significantly to the total spectrum. This suggests 
that modest changes to the outflow kinematics might boost 
the wind continuum and produce emergent spectra with 
weak or absent Balmer absorption edges. 
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Figure 7. Synthetic optical spectra from model A computed for 
sightlines of 10, 27.5, 45, 62.5 and 80 degrees. In these plots the 
flux is divided by a polynomial fit to the underlying continuum 
redward of the Balmer edge, so that line-to-continuum ratios and 
the true depth of the Balmer jump can be shown. 



Wavelength (A) 

Figure 8. Total packet-binned spectra across all viewing angles, 
in units of monochromatic luminosity. The thick black line shows 
the total integrated escaping spectrum, while the green line shows 
disk photons which escape without being reprocessed by the wind. 
The red line show the contributions from reprocessed photons. 
Recombination continuum emission blueward of the Balmer edge 
is already prominent relative to other wind continuum processes, 
but is not sufficient to fill in the Balmer jump in this specific 
model 
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5 A REVISED MODEL OPTIMIZED FOR 

OPTICAL WAVELENGTHS 

The benchmark model discussed in section |4] was originally 
designed to reproduce the wind-formed lines seen in the UV 
spectra of high-state CVs. As we have seen, this model does 
produce some observable optical emission. We can now at¬ 
tempt to construct a model that more closely matches the 
observed optical spectra of CVs. 

Specifically, we aim to assess whether a revised model 

can: 

• account for all of the lines we see in optical spectra of 
CVs while preserving the UV behaviour; 

• produce single-peaked Balmer emission lines; 

• generate enough of a wind-formed recombination con¬ 
tinuum to completely fill in the disk’s Balmer absorption 
edge for reasonable outflow parameters. 

The emission measure of a plasma is directly propor¬ 
tional to its density. The simplest way to simultaneously 
affect the density in the wind (for fixed mass loss rate), as 
well as the velocity gradients, is by modifying the poloidal 
velocity law. Therefore, we focus on just two kinematic vari¬ 
ables f section |3.1l) : 

• the acceleration length, , which controls the distance 
over which the wind accelerates to i Vao', 

• the acceleration exponent, a, which controls the rate at 
which the poloidal velocity changes near 

The general behaviour we might expect is that outflows 
with denser regions near the wind base - i.e. winds with 
larger and/or larger a - will produce stronger optical 
emission signatures. However, this behaviour may be mod¬ 
erated by the effect of the increasing optical depth through 
this region, which can also affect the line profile shapes. In 
addition, modifying also increases the emission volume. 
Based on a preliminary exploration of models with different 
kinematics, we adopt the parameters listed in table[T]for our 
‘optically optimized’ model (model B). 

5.1 Synthetic Spectra 

Fig.|9]shows the UV and optical spectra for the optically op¬ 
timized model for the full range of inclinations. As expected, 
the trend from absorption to emission in the optical is again 
present, but in this revised model we produce emission lines 
in the entire Balmer series at high inclinations, as well as 
the observed lines in He ll and He l. This can be seen more 
clearly in the continuum-normalized spectrum in Fig. [10] 

Two other features are worth noting in the optical spec¬ 
trum. First, the collisionally excited Ca ll emission line at 
3934 A becomes quite prominent in our densest models. Sec¬ 
ond, our model predicts a detectable He ll recombination line 
at 3202 A. This is the He equivalent of Paschen [3 and should 
be expected in all systems that feature a strong He ll A4686 
line (the He equivalent of Paschen a). This line is somewhat 
unfamiliar observationally, because it lies bluewards of the 
atmospheric cut-off, but also redwards of most ultraviolet 
spectra. 

Our models do not exhibit P-Cygni profiles in the op¬ 
tical lines. This is perhaps not surprising. LK02 and SV93 
originally designed such models to reproduce the UV line 


profiles. Thus, most of the wind has an ionization param¬ 
eter of log!/ 2 (see Fig. [5|). This means H and He are 
fully ionized throughout much of the wind and are success¬ 
ful in producing recombination features. However, the line 
opacity throughout the wind is too low to produce notice¬ 
able blue shifted absorption. We suspect that the systems 
that exhibit such profiles must possess a higher degree of 
ionization stratification, although the lack of contemporary 
observations means it is not known for certain if the P-Cygni 
profiles in UV resonance lines and optical H and He lines ex¬ 
ist simultaneously. Ionization stratification could be caused 
by a clumpy flow, in which the ionization state changes due 
to small scale density fluctuations, or a stratification in den¬ 
sity and ionizing radiation field over larger scales. Invok¬ 
ing dumpiness in these outflows is not an unreasonable hy- 
poth esis. Theories of line-driven winds predict an unstabl e 
flow llMacGregor et all 19791 : lOwocki fc Rvbickil984ll985l i. 
and simulation s of CV disk winds also produce density in¬ 
homogeneities (jProga et al.lll99^ . l2002h . Tentative evidence 
for clumping being dire ctly related to P-Cy gni optical lines 
comes from the fact that jPrinia et al] (l200Cll l found the dwarf 
nova BZ Cam’s outflow to be unsteady and highly mass- 
loaded in outburst, based on observations of the UV reso¬ 
nance lines. This system has also exhibited P-Cygni profiles 
in He l A5876 and Ha when in a high-state ijPatterson et al.l 
1 19961 : iRingwald fc Navloi]|l99^ . The degree of ionization 
and density variation and subsequent line opacities may be 
affected by our model parameters and the specific parame- 
terisation we have adopted. 

In the UV, the model still produces all the observed 
lines, and deep P-Cygni profiles are produced in the normal 
resonance lines, as discussed in section 4.2. However, the UV 
spectra also display what is perhaps the biggest problem 
with this revised model, namely the strength of resonance 
line emission at low and intermediate inclinations. In order 
to generate strong optical wind signatures, we have adopted 
wind parameters that lead to very high densities at the base 
of the wind (ue ~ 10^^ — 10^"* cm“®). This produces the 
desired optical recombination emission, but also increases 
the role of collisional excitation in the formation of the UV 
resonance lines. This explains the pronounced increase in 
the emission component of the CiV A1550 resonance line, for 
example, relative to what was seen in the benchmark model 
(compare Figures [6] and [9||. The strength of this component 
in the revised model is probably somewhat too high to be 
consistent with UV observations of high-state CVs (see e.g. 
Long et al. 1991, 1994; Noebauer et al. 2010). 


5.2 Continuum Shape and the Balmer Jump 

The wind now also has a clear effect on the continuum 
shape, as shown by Fig. 1111 In fact, the majority of the 
escaping spectrum has been reprocessed in some way by the 
wind, either by electron scattering (the wind is now mod¬ 
erately Thomson-thick), or by bound-free processes. This is 
demonstrated by the flatter spectral shape and the slight 
He photoabsorption edge present in the optical spectrum 
(marked in Fig. 1101) . This reprocessing is also responsible 
for the change in continuum level between models A and B. 
In addition. Figures 1101 and 1111 clearly demonstrate that 
the wind produces a recombination continuum sufficient to 
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Figure 9. UV (left) and optical (right) synthetic spectra for model B computed at sightlines of 10, 27.5, 45, 62.5 and 80 degrees. 
Model A is shown in grey for comparison. The inset plots show zoomed-in line profiles for He ii A1640 and Ho. The Balmer and He are 
double-peaked, albeit with narrower profiles. Strong He ii A4686 emission can be seen, as well as a trend of a deeper Balmer jump with 
decreasing inclination. 
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completely fill in the Balmer jump at high inclinations0 This 
might suggest that Balmer continuum emission from a wind 
can be important in shaping the Balmer jump region, as 
originally suggested by Knigge et al. (1998b; see also Has- 
sall et al. 1985). 

It should be acknowledged, however, that the Balmer 
jump in high-state CVs would naturally weaken at high in- 
clinations due to limb darkening effects (iLa Douslll989bl l^. 
Although we include a simple limb darkening law which af¬ 
fects the emergent flux at each inclination, we do not include 
it as a frequency dependent opacity in our model. As a re¬ 
sult, the efficiency of filling in the Balmer jump should really 
be judged at low and medium inclinations, where, although 
prominent, the recombination continuum does not overcome 
the disk atmosphere absorption. In addition, this effect could 
mean that any model which successfully fills in the jump at 
low inclinations could lead to a Balmer jump in emission at 
high inclinations. In any case, to properly understand this 
phenomenon, a fully self-consistent radiative transfer calcu¬ 
lation of both the disk atmosphere and connected wind is 
required. 


5.3 Line Profile Shapes: Producing Single-Peaked 
Emission 

Fig. [12] shows how the Ho profile changes with the kine¬ 
matics of the wind for an inclination of 80°. The main pre¬ 
diction is that dense, slowly accelerating wind models pro¬ 
duce narrower emission lines. This is not due to radial ve¬ 
locity shear. As stated by MC96, that mechanism can only 
work if poloidal and rotational velocity gradients satisfy 
{dvi/dr)/{dv^/dr) > 1; in our models, this ratio is always 
< 0.1. Instead, the narrow lines predicted by our denser wind 
models can be traced to the base of the outflow becoming 
optically thick in the continuum, such that the line emission 
from the base of the wind cannot escape to the observer. 
In such models, the ‘line photosphere’ (the r ~ 1 surface 
of the line-forming region) moves outwards, towards larger 
vertical and cylindrical distances. This reduces the predicted 
line widths, since the rotational velocities - which normally 
provide the main line broadening mechanism at high incli¬ 
nation - drop off as 1/r. This is not to say that the MC96 
mechanism could not be at work in CV winds. For example, 
it would be worth investigating alternative prescriptions for 
the wind velocity field, as well as the possibility that the out¬ 
flows may be clumped. An inhomogeneous flow (which has 
been predicted in CVs; see section 5.2) might allow large ra¬ 
dial velocity shears to exist while still maintaining the high 
densities needed to produce the required level of emission. 
However, such an investigation is beyond the scope of the 
present paper. 

In our models, single-peaked line profiles are produced 
once the line forming region has been pushed up to ~ 
10^^ cm (~ 150 Rwd) above the disk plane. This num¬ 
ber may seem unrealistically large, but the vertical extent 


^ Note that the apparent absorption feature just redward of the 
Balmer jump in these models is artificial. It is caused by resid¬ 
ual line blanketing in the stellar atmospheres, which our models 
cannot fill in since they employ a 20-level H atom. 



Wavelength (A) 

Figure 10. Synthetic optical spectra from model B computed 
for sightlines of 10, 27.5, 45, 62.5 and 80 degrees. Model A is 
shown in grey for comparison. In these plots the flux is divided 
by a polynomial fit to the underlying continuum redward of the 
Balmer edge, so that line-to-continuum ratios and the true depth 
of the Balmer jump can be shown. 

of the emission region is actually not well constrained ob- 
servationally. In fact, multiple observations of eclipsing NLs 
show that the Ha line is only moderately eclipsed compared 
to the continuum (e.g. Baptista et al. 2000; Groot et al. 
2004; see also section 5.4), implying a significant vertical ex¬ 
tent for the line-forming region. This type of model should 
therefore not be ruled out a priori, but this specific model 
was not adopted as our optically optimized model due to its 
unrealistically high continuum level in eclipse. 

5.4 Sensitivity to Parameters 

This revised model demonstrates that one can achieve a 
more realistic optical spectrum by altering just two kine¬ 
matic parameters. However, it may also be possible to 
achieve this by modifying other free parameters such as 
the opening angles of the wind and the inner and 
outer launch radii. For example, increasing the mass loss 
rate of the wind increases the amount of recombination emis¬ 
sion (which scales as p^), as well as lowering the ionization 
parameter and increasing optical depth through the wind. 
Wider launch radii and opening angles lead to a larger emit¬ 
ting volume, but this is moderated by a decrease in density 
for a fixed mass-loss rate. We also note that the inner ra¬ 
dius of 4 Rwd adopted by SV93 affects the emergent UV 
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Wavelength (A) 

Figure 11. Total packet-binned spectra across all viewing angles, 
in units of monochromatic luminosity. The thick black line shows 
the total integrated escaping spectrum, while the green line shows 
disk photons which escape without being reprocessed by the wind. 
The red line show the contributions from reprocessed photons. 
In this denser model the reprocessed contribution is significant 
compared to the escaping disk spectrum. The Balmer continuum 
emission is prominent, and the wind has a clear effect on the 
overall spectral shape. 



Figure 12. Ha line profiles, normalized to 1, plotted in velocity 
space for three models with varying kinematic properties, com¬ 
puted at an inclination of 80°. The benchmark model and the 
improved optical model described in section 6 are labeled as A 
and B respectively, and a third model (X) which has an increased 
acceleration length of Ry = 283.8 Rwd, s,nd a = 4 is also shown. 
The a:-axis limits correspond to the Keplerian velocity at 4Rw d , 
the inner edge of the wind. We observe a narrowing of the lines, 
and a single-peaked line in model X. This is not due to radial 
velocity shear (see section 5.3). 


Spectrum seen at inclinations < Omin as the inner disk is 
uncovered. This causes less absorption in the UV resonance 
lines, but the effect on the optical spectrum is negligible. We 
have verified this general behaviour, but we suggest that fu¬ 
ture work should investigate the effect of these parameters 
in more detail as well as incorporating a treatment of clump¬ 


ing. If a wind really does produce the line and continuum 
emission seen in optical spectra of high-state CVs, then un¬ 
derstanding the true mass loss rate and geometry of the 
outflow is clearly important. 


5.5 Comparison to RW Tri 

Fig. shows a comparison of the predicted out-of-eclipse 
and mid-eclipse spectra against observations of the high- 
inclination nova-like RW Tri. The inclination of RW Tri is 
some what uncertain, with est imates including 70.5° ( SmakI 
Il995h . 7 5° dCroot et al.ll2004h . 80° (|Longmore et al.l 1981 ) 
and 82° ([Frank &: Kindll98lh . Here, we adopt i = 80°, but 
our qualitative conclusions are not particularly sensitive to 
this choice. We follow LK02 is setting the value of rdisk 
(the maximum radius of the accretion disk) to 34.3 Rwd- 
When compared to the semi-major axis of RW Tri, this 
value is perhaps lower than one m ight typically expect for 
NLs llHarrop-Allin fc Warner 1996ll. However, it is consis¬ 
tent with values inferred by Rutten et al.l (Il992ll . We em¬ 
phasize that this model is in no sense a fit to this - or any 
other - data set. 

The similarity between the synthetic and observed 
spectra is striking. In particular, the revised model pro¬ 
duces strong emission in all the Balmer lines, with line- 
to-continuum ratios comparable to those seen in RW Tri. 
Moreover, the line-to-continuum contrast increases during 
eclipse, as expected for emission produced in a disk wind. 
This trend is in line with the observations of RW Tri, and it 
has also been seen in other NLs, incl uding members of the 
SW Sex class (iNeustroev et al]l201ll 'l. As noted in section 
5.2, the majority of the escaping radiation has been repro¬ 
cessed by the wind in some way (particularly the eclipsed 
light). 

However, there are also interesting differences between 
the revised model and the RW Tri data set. For example, 
the model exhibits considerably stronger He ll features than 
the observations, which suggests that the overall ionization 
state of the model is somewhat too high. As discussed in 
section 5.3, the optical lines are narrow, but double-peaked. 
This is in contrast to what is generally seen in observations 
of NLs, although the relatively low resolution of the RW 
Tri spectrum makes a specific comparison difficult. In order 
to demonstrate the double-peaked nature of the narrower 
lines, we choose not to smooth the synthesized data to the 
resolution of the RW Tri dataset. If the data was smoothed, 
the Ha line would appear single-peaked. 


6 CONCLUSIONS 

We have investigated whether a disk wind model designed 
to reproduce the UV spectra of high-state CVs would also 
have a significant effect on the optical spectra of these sys¬ 
tems. We find that this is indeed the case. In particular, the 
model wind produces H and He recombination lines, as well 
as a recombination continuum blueward of the Balmer edge. 
We do not produce P-Cygni profiles in the optical H and He 
lines, which are seen in a small fraction of CV optical spec¬ 
tra. Possible reasons for this are briefly discussed in section 
5.2. 

We have also constructed a revised benchmark model 
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Wavelength (A) 



Wavelength (A) 

Figure 13. Top Panel: In and out of eclipse spectra of the high inclination NL RW Tri. Bottom Panel: In and out of eclipse synthetic 
spectra from model B. The artificial ‘absorption’ feature just redward of the Balmer jump is caused for the reasons described in section 
5.2. 


which is designed to more closely match the optical spectra 
of high-state CVs. This optically optimized model produces 
all the prominent optical lines in and out of eclipse, and 
achieves reasonable verisimilitude with the observed opti¬ 
cal spectra of RW Tri. However, this model also has signif¬ 
icant shortcomings. In particular, it predicts stronger-than- 
observed He ll lines in the optical region and too much of a 
collisionally excited contribution to the UV resonance lines. 

Based on this, we argue that recombination emission 
from outflows with sufficiently high densities and/or opti¬ 
cal depths might produce the optical lines observed in CVs, 
and may also fill in the Balmer absorption edge in the spec¬ 
trum of the accretion disk, thus accounting for the absence 
of a strong edge in observed CV spectra. In section 5.3, we 
demonstrate that although the double peaked lines narrow 
and single-peaked emission can be formed in our densest 
models, this is not due to the radial velocity shear mech¬ 
anism proposed by MC96. We suggest that ‘clumpy’ line- 
driven winds or a different wind parameterization may nev¬ 
ertheless allow this mechanism to work. We also note the 
possibility that, as in our denser models, the single-peaked 
lines are formed well above the disk, where rotational veloc¬ 
ities are lower. 

It is not yet clear whether a wind model such as this 
can explain all of the observed optical features of high-state 
CVs - further effort is required on both the observational 
and modelling fronts. However, our work demonstrates that 
disk winds matter. They are not just responsible for creating 
the blue-shifted absorption and P-Cygni profiles seen in the 
UV resonance lines of high-state CVs, but can also have a 
strong effect on the optical appearance of these systems. In 


fact, most of the optical features characteristic of CVs are 
likely to be affected - and possibly even dominated - by 
their disk winds. Given that optical spectroscopy plays the 
central role in observational studies of CVs, it is critical to 
know where and how these spectra are actually formed. We 
believe it is high time for a renewed effort to understand 
the formation of spectra in accretion disks and associated 
outflows. 
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